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The Hong-Ou-Mandel effect is generahzed to a configuration of n bosons prepared in the n input 
ports of a BeU muhiport beam splitter. We derive a strict suppression law for most possible output 
events, consistent with a generic bosonic behavior after suitable coarse graining. 
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The Hong-Ou-Mandel (HOM) effect [J is an impres- 
sive manifestation of the bosonic quantum nature of pho- 
tons. In the original experiment, two identical photons 
are sent simultaneously (within their coherence time) 
through the two input ports of a balanced beam splitter. 
Since no interaction between the photons takes place, 
one would intuitively expect the photons to propagate 
independently and not presume any correlations in the 
number of photons measured at both output ports. Sur- 
prisingly, the photons always leave the setup together, 
but never exit at different ports. Such coincident events 
at both output ports are completely suppressed. 

This effect is used in many applications: The visibility 
of the dip in the coincident detection rate provides a char- 
acterization for the indistinguishability of two photons 
[21 [3], and therewith for the quality of photon sources. 
HOM setups are used to project photons onto the max- 
imally entangled l^*") Bell-state, and consequently to, 
both, create and detect such states [1]. This is used 
for example in entanglement swapping protocols 5J and 
quantum metrology 6 . Furthermore, the nondeterminis- 
tic gate operations in linear optics quantum computation 
[7] are based on the HOM effect. 

It is suggestive to generalize the HOM setup for more 
than two photons and more than two input or output 
ports. Indeed, the enhancement of events with all par- 
ticles in one port - bunching events - has been observed 
experimentally when several photons enter each of the 
two modes of an unbiased (i.e., balanced) two-port beam 
splitter [5J 12] . For a specially designed biased setup of 
three particles and three input ports, the suppression 
of coincident events was shown [10]. In the case of a 
Bell multiport beam splitter [TTJ [T^] which redistributes 
n incoming particles to n ports in an unbiased way, it is 
known that coincident events are suppressed when n is 
even [T5] . 

All these results imply important applications, from 
the creation and detection of multipartite qudit- 
entangled states [HI [15] , over the implementation of en- 
tanglement swapping protocols for many particles and 



the design of efficient quantum gates for qudits [TB], to 
the experimentally controlled transition from indistin- 
guishability to distinguishability for many identical par- 
ticles However, we still lack a comprehensive under- 
standing of the n-particle, n-port generalization of the 
HOM effect, since the complexity of such a scattering 
problem scales very unfavorably with n: The number of 
interfering amplitudes as well as that of possible output 
events grow faster than exponentially. Hence, a detailed 
analysis of individual output events is prohibitive, and 
needs to be substituted by statistical considerations. 

This is the purpose of the present Letter, where we 
present a general study of the probabilities of all pos- 
sible output events of the Bell multiport beam splitter. 
Our treatment enables a general understanding of mul- 
tiparticle interference effects, as well as on the average 
behavior of bosons. It hence unifies previous experimen- 
tal and theoretical work on multiport beam splitters, and 
opens up new perspectives for the experimental verifica- 
tion and exploitation of bosonic multiparticle behavior. 

In the following we denote arrangements of n parti- 
cles in the n modes by a vector s — (si, S2, • • • Sn), with 
Sfc the number of particles in the output mode k, and 
X]"=i ~ For distinguishable, noninterfering parti- 
cles the probability for a certain arrangement s reads 

1 n\ 

We call this situation "classical", since, due to the lack 
of interference between the particles, probabilities are 
summed instead of amplitudes, and simple combinatorics 
applies. Hence, coincident events, i.e. Sc = (1,1,... 1), 
are realized with probability n\ln^. Bunching events, 
with all particles at one output mode k, correspond to 
Sk — n and thus to Sb = (0, 0, .., n, ..0). They are realized 
with probability 1/n" and, hence suppressed by a factor 
of n! with respect to the coincident events. For large n, 
both events are highly unlikely, extreme cases. 

The analogous problem with identical quantum par- 
ticles is best formulated in second quantization. Since 



2 



applications of our study are feasible with today's opti- 
cal technologies [IT], we focus here on bosons with the 
following commutation relations for the creation and an- 
nihilation operators for the respective ports: 



0. 



(2) 



The initial state reads j^*) ~ Y[i=i'^l\^)- Input port 
creation operators a| are mapped to output creation 
operators via a unitary matrix U [il8 , such that 

= X]fc=i UjkO'\- Formally, the unbiased Bell multiport 
beam splitter under consideration corresponds to the uni- 
tary operation given by the Fourier matrix, defined for 
any dimension n by Ujk = e~^^^~^'''^^~^^ / \/n. 

The possible states with fixed particle number per port 
after the scattering process read 
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|0). 



(3) 



In order to describe the event probability of a given ar- 
rangement s, we define a vector d of length n with entries 
that specify each particle's output port. It is constructed 
by concatenating Sj times the port number j: 



(4) 



e.g., for the arrangement s — (2, 1, 0, 2, 0), we find d{s) = 
(1, 1, 2, 4, 4). Therewith we can write the transition prob- 
ability to a specific output arrangement s 



cr£P„ 3 = 1 



,(5) 



where P„ denotes the set of all permutations of {1, .., 7i}. 
For coincident output states, we have dj{sc) = j; i.e., 
the overlap ([5| becomes the permanent of the matrix U 
[ini [2D] . The best known algorithm to compute the am- 
plitude ([5| scales exponentially in n \21\ . In information- 
theoretic terms, the evaluation hence remains an NP- 
complete problem [T3], despite the symmetry of U. 

Notwithstanding the apparent complexity of the prob- 
lem, it is possible to exploit the symmetry of the ma- 
trix U to formulate a powerful law (with its proof given 
in the appendix) : Events characterized by s are strictly 
suppressed if the sum of vector-components d;(s) is not 
dividable by n: 

n 

Q(s-) :=Mod(^d,(s),n) ^0 ^ (*|$(.s)) = 0. (6) 
1=1 

Consider, e.g., n = 6 and si = (2,1,2,1,0,0): One 
immediately finds Q{si) — 2, and this event is hence 
strictly suppressed. Unexpectedly though, the event 
5*2 = (0, 1, 2, 0, 2, 1), which is obtained from si by simple 



permutation, gives Q{s2) = 0, and is actually enhanced 
by a factor larger than 7 as compared to the classical 
event probability (also see Table |ll]). Note that the eval- 
uation of ([6]) scales linearly with n, and that - as shown 
below - our suppression law applies for most output ar- 
rangements. It thus largely characterizes the general sta- 
tistical behavior of the Bell multiport beam splitter, in a 
easily evaluable manner. 

For a more detailed insight in the predictions of Eq. (|6| , 
we first have to identify classes of final states that occur 
with equal probability. In the classical case, the realiza- 
tion probability of any arrangement s remains invariant 
under permutation of the output ports s^. Hence we 
can define classical equivalence classes which identify ar- 
rangements related to each other by permutation. The 
amplitude ([s]), however, is not invariant under arbitrary 
permutations of the Sk] i.e., two classically equivalent 
arrangements are not necessarily quantum mechanically 
equivalent. Only cyclic and anticyclic permutations leave 
Eqs. (5[6) invariant. This allows us to define a quantum 
equivalence relation between arrangements, and the as- 
sociated quantum equivalence classes. To estimate the 
number of suppressed arrangements that are predicted 
by ([6]), let us assume that the Q{s) are uniformly dis- 
tributed in the interval [0, . . . , n — 1] for the ensemble of 
events s. Then the probability to find a suppressed ar- 
rangement is given by the weight of nonvanishing values 
of Q{s), i.e., by (1 — n)/n = 1 — 1/n. The number of 
equivalence classes and the number of suppressed out- 
put arrangements, shown in Table |Tj for n = 2. .14, also 
confirm this scaling numerically. 

The unsuppressed arrangements are listed in Table |llj 
for n — 2. .6, together with their quantum enhancement, 
i.e., the ratio of quantum-to-classical event probability. 
The suppression of coincident events for even n derived 
in [13 is covered by Eq. ([6]) as a special case: For such 
events, d{sc) = (1, 2, 3, .., n); hence dj(s) = n(n-|-l)/2 
is never dividable by n for even n. 

The implications of Eq. (|6| are rather counterintu- 
itive. On the one hand, bunching events are enhanced 
by a factor of nl with respect to the classical case, as ex- 
pected due to the bosonic nature of particles: This favors 
states with many particles in few occupied ports. On the 
other hand, the number of particles in one port, or the 
number of occupied ports is not a direct indicator for 
the enhancement or the suppression of a certain event. 
For example, one intuitively expects events of the type 
s = {n — 1, 1, 0, .., 0) to be enhanced due to the bosonic 
nature of the particles, while they actually turn out to 
be strictly suppressed. Thus, at the level of the event 
probabilities of single arrangements, interference effects 
dominate, and the bosonic nature of the particles is not 
apparent at all. 

It is, however, possible to recover a general bosonic 
behavior by grouping many final arrangements in larger 
classes which are characterized, e.g., by the number of 
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TABLE I: Number of classical equiva- 
lence classes (A^ciass), quantum equiv- 
alence classes (A'quantum), classes that 
therewithin fulfill the law Q (Maw), 
and suppressed classes which are not 
predicted by Eq. |6] (Mupp). 



n s 


Enhancement 


1 


3 (003) 


6 




(111) 


3/2 




4 (0004) 


24 




('0202") ('0121") 


8/9 




5 (00005) 


120 


c3 

Oh 


(00131), (01103) 


15/2 


o 


(00212), (01022) 


10/3 




(11111) 


5/24 


|10 


6 ('000006'! 


720 




(002004), (000141), 






(010104) , (000303), 


1 144/5 


14 


(001032), (000222) 






(020202), (001113), (012021) 


36/5 





TABLE IL Nonsuppressed output 
states, together with the correspond- 
ing quantum enhancement, i.e., the ra- 
tio of quantum to classical event prob- 
ability. 
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TABLE III: (color online) Quantum 
enhancement as a function of the num- 
ber of particles (vertical axis), and of 
the number of occupied ports (hori- 
zontal axis). 



occupied ports k, or by the number m of particles in one 
port. The event probability for such a class is given by 
the sum of the probabilities of the single events that per- 
tain to the class. Very generally, one expects that, for 
bosons, quantum states with large occupation numbers 
are favored. This general behavior is also reflected by our 
formalism: According to (|5|, the probabilities Pqm{s) are 
given in terms of a sum over permutations of scattering 
amplitudes, i.e., over complex numbers of equal modulus 
(products of matrix-elements of Ujk)- Since these num- 
bers typically have different phases, they tend to add up 
destructively. However, all Sj\ permutations a that in- 
terchange the Sj particles that exit in port j leave the 
scattering amplitudes invariant, so that Sj\ terms in the 
sum have equal phases and add up constructively. This 
motivates the following approximation for the transition 
probability ([s]): 
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We show the probability distribution for the number of 
occupied ports, for the classical calculation ([T]), for the 
bosonic quantum case ( [5|), and for our approximation 
([7]), for n ^ 14, in Figure 111 The distributions' expecta- 
tion values correspond to the average number of occupied 
ports. This value can be shown to grow linearly with the 
number of particles n. However, the growth rate in the 
bosonic case is approx. 0.50, smaller than the value of 
approx. 0.63 found in the classical case. As expected, 
bosons always tend to occupy less output ports than in 
the classical case, for any n. Furthermore, Figure [T] shows 
that the approximation ([t]) predicts the actual outcome 



very well for most k, and only fails for events with almost 
all or almost no sites occupied. This is easily understood, 
since, for very small or very large fc, few distinct equiv- 
alence classes contribute to these event groups. Then, 
again interference dominates the event probability, rather 
than bosonic behavior. Therefore, a few general suppres- 
sion effects that follow from ^ persist at this level, for 
extreme values of k: Events with n — 1 occupied ports 
are suppressed for odd n, coincident events are forbid- 
den for even n. If n is prime, there will be never exactly 
two occupied ports i,j: The sum in ([6]) then becomes 
k ■ i + in — k)j with k the number of particles in port i. 
The result is never dividable by n for < fc < n. These 
effects and the quantum enhancement of the classes with 



k occupied ports are visualized in Table HI 



Also the event probability for a given number of par- 
ticles in one single port is well described by our esti- 
mate ([t]). For 14 particles, the probability distribution is 
shown in Figure [2j Again, we see a dramatic difference 
between the classical and quantum case, especially for 
the probability to find a large number of particles in one 
port. 

We have generalized the HOM effect to n particles and 
n ports on two different levels: Interference effects inhibit 
the realization of most possible events for single transi- 
tion amplitudes, while general statistical characteristics 
with smooth bosonic behavior emerge that are efficiently 
approximated by Eq. ([T]). On the fine as well as on 
the coarse grained scale, however, quantum and classical 
transmission probabilities differ dramatically. 

In order to verify our above theoretical findings, and to 
statistically characterize the indistinguishability of many 
photons, single-photon counting detectors are required 
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FIG. 1: (color online) Event probability for a given number of 
occupied ports, for n = 14. (Red) rectangles denote classical 
combinatorics, (blue) triangles the quantum mechanical prob- 
ability distribution, and (black) circles our (bosonic) estimate 
for the quantum result. The inset shows the same distribution 
on a linear scale. Note that events with 14 occupied ports are 
strictly suppressed in the quantum case. 
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Particles per port 

FIG. 2: (color online) Probability to find exactly k particles 
(horizontal axis) in one port, for n — 14. (Red) rectangles 
denote the classical, (blue) triangles the quantum calculation, 
and (black) circles the estimate ([7|. The inset shows the 
distribution on a linear scale for small k. Note that events 
with 13 particles in one port are totally suppressed in the 
quantum case. 



in the experiment. Since such detectors are not stan- 
dard equipment yet, let us stress that also a more coarse 
grained measurement with bucket detectors which do not 
count the number of simultaneously arriving photons ex- 
hibits a strong quantum signature in the event statistics, 
as clearly spelled out in Figure [T] 
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Appendix. Each of the n\ terms in the sum in Eq. 
can be written as an n-th root of unity. Hence, Eq. 
turns into X]fc=o '^ke^^'^, where the Ck are natural 
numbers which give the cardinality of the following sets. 



.(s) = ia|e„,,» = ^d(.s)(a(0-rmodnL(8) 



1=1 



with Cr — \ur{s)\. The sum corresponds to the position of 
the barycenter of the set of points {cke^~^'^\k G {1, .., n}} 
in the complex plane. We set Q = mod (^"^^^ d/(s), n), 
and define an operation 7 which acts on permutations 
such that 7((7)(fc) = cr(A:) + 1 mod n. It is immediate 
that en,sili<y)) = 0™,s-(o-) + Q- Thus, if Q 7^ 0, the 
repeated application of 7 gives us a bijection between all 
pairs of Wj+a Q, for a G {0, 1, ..,n — 1}. Hence, we find 



Vr e {0, 



l},Va e N : Cr+a Q 



(9) 



Therefore, if Q 7^ 0, the set of points {cfee*^'^|fc g 
{l,..,n}} describes several interlaced polygons centered 
at the origin, ensuring that the sum vanishes. On the 
other hand, it is also possible for the barycenter of the 
structure spanned by the to lie in the origin, even 
though the set of points is not described by polygons. 
Therefore the reverse of the law ^ does not hold. 
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